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c*| , Abstract 

An empirical principle for the construction of a linear relationship between 
^ ■ the total angular momentum and squared-mass of baryons is proposed. In or- 

der to examine linearity of the trajectories, a rigorous least-squares regression 
analysis was performed. Unlike the standard Regge-Chew-Frautschi approach, 
the constructed trajectories do not have non-linear behaviour. A similar reg- 
ularity may exist for lowest-mass mesons. The linear baryonic trajectories 
are well described by a semi-classical picture based on a spinning relativistic 
string with tension. The obtained numerical solution of this model was used 
to extract the (di)quark masses. 
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1 Introduction 



It has been accepted for a long time that hadrons from the same family lie on Regge 
trajectories (the so-called Chew-Frautschi conjecture [1]), i.e. 

J = a(0)+aM 2 , (1) 

where J is the total angular momentum and M is the mass of a hadron. The 
intercept a(0) depends on hadron type, but the slope a is approximately the same 
for all hadrons. Such a relationship between J and M 2 , also known as the principle of 
exchange degeneracy, is usually interpreted as a manifestation of the linear potential 
of the strong forces between constituent quarks. Recently this picture was discussed 
in terms of a relativistic diquark model based on a spinning string with a constant 
tension [2-4] (see also references in [5]). 

Experimentally, the validation of Eq. ((TJ remains to be a difficult problem since 
the experimental data are scarce. According to a recent classification of hadrons on 
the (J, M 2 ) plane [6], the overwhelming majority of the trajectories are supported 
by a few data points only. For mesons, there exist five trajectories with three data 
points, while other trajectories were hypothesized from two or one data points. For 
baryons, only one trajectory was constructed from four data points and five trajec- 
tories were supported by three data points. Other trajectories were hypothesized 
from the study of either one or two data points. Recently, it has been noted that 
only 9% (14%) of all trajectories in the mesonic (baryonic) sector are linear [7]. A 
non-linear character of the Regge trajectories has also been pointed out in [4,8,9]. 

In this article we are not going to scrutinize the linearity of the Regge trajectories; 
it is rather clear that the vast majority of such trajectories are indeed non-linear 
and there is no necessity to re-analyse this fact again. Instead, using the most 
recent PDG data [10], we would like to note that it is possible to find a prescription 
which could allow the construction of the trajectories which can be classified as 
being perfectly linear and span over a significant number of known baryons. In 
order to illustrate this, we use a rigorous weighted least-squares regression, which 
is often missing in theoretical papers on this subject. In Sect. 121 we will discuss 
our principles for the construction of a linear relationship between the hadronic 
mass squared M 2 and the total angular momentum J. In Sect. our empirical 
observation is discussed in more details. In particular, we will explain what would 
happen if the requirements proposed in Sect. |2] are removed or weakened. In Sect. 0] 
and appendices, we will attempt to use our approach for validation of a relativistic 
diquark model based on a spinning string, which is often discussed in the literature 
(see [4] for a recent discussion). Finally, we will briefly discuss the mesonic sector 
in Sect. |S] 
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2 Linear trajectories for baryons 



As a starting point in the construction of baryonic trajectories, we will consider: (a) 
only stable or strongly decaying particles from a certain family; (b) only hadrons 
of the same family with a smallest mass, M m i n j, for a given fixed J and space 
parity P; (c) in addition to (a) and (b), we will use only such baryons if there are 
no other baryons with smaller masses and with the opposite space parity. At this 
moment, we will use the assumptions above without any theoretical justification. 
Below we will show that, without any exception, such requirements are sufficient for 
the construction of perfectly linear baryonic trajectories. 

As a leading principle for the selection of hadrons, we will ignore particles with 
poor evidence of existence i.e. with one star in accordance with the PDG classifi- 
cation [10]. Furthermore, in several cases when PDG quotes only mass ranges, we 
will use the central values for the masses and their experimental errors as they are 
given by the most accurate and recent measurements 1 . 

We will attempt to describe the mass dependence on J using a linear parame- 
terisation similar to Eq. (|T|): 

MLn,j=Po+PiJ, (2) 

where p$ and p\ are constants. Hereafter we will use the notation M for M min j, 
unless otherwise stated. For convenience of a linear-regression analysis, we prefer 
to express M 2 as a function of J since the angular momentum does not have ex- 
perimental uncertainties. The linear least-squares regression analysis is described in 
Appendix A. 

TV baryons. First, let us consider N baryons. Figure Q shows all PDG N- 
baryons on the (J,M 2 ) plane (the filled and open symbols). There is a clear linear 
trend for lowest-mass baryons with P = +1 shown by the filled circles (the PDG 
names of such baryons are indicated). Since such baryons fully satisfy the criteria 
proposed at the beginning of this section, we will use them for a linear least-squares 
regression. The result of this regression is shown with the solid line. The trajectory 
is remarkably straight: the small value of x 2 /ndf shown in Fig. [T] fully supports the 
linear fit. 

To check a possible non-linearity of the iV-baryon trajectory, a fit was performed 
using a second-order polynomial function. The inclusion of additional term to the 
fitting function leads to a small value (~ 0.005) of the parameter responsible for the 
quadratic term. However, the quality of such fit characterized by x 2 /ndf = 0.15/1 
does not significantly improve. Two dashed lines shown in Fig. ^ indicate a 95% 
confidence-level region for the linear regression, assuming that the experimental 
uncertainties on the masses are normally distributed. This illustrates the reliability 
of the linear relationship between the M 2 and J values. From the two tests above, 

1 This can be found in the computer files located at the PDG web page [11]. 
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it can be concluded that there is no evidence for non-linearity of the iV-baryon 
trajectory 

It should be noted that we do not include the iV(1990)(-F17)** baryon at J = 7/2 
to the fit. In principle, the quality of the fit will not suffer if the quoted PDG mass 
(~ 1990 MeV) is used. However, if one uses the most recent measurement given by 
the PDG [11], then the iV(1990)(F17) point will move up to the location indicated by 
the small symbol at 2086 ±28 MeV (Fig. [1}. In this case, this baryon will overlap on 
the (J, M 2 ) plane with a better established P = —1 baryon iV(2190)(Gl7), which has 
a mass of 2190^ MeV (again in accordance with the latest measurement). In this 
case, iV(1990)(F17) should not be used for the fit due to the "lowest-mass" exclusion 
principle. Since there is no any objective criteria for inclusion (or exclusion) of non- 
well established iV(1990)(F17) baryon to (from) the fit, it was decided to exclude 
it. 

The situation with the lowest-mass P = —1 N baryons, shown in Fig. ^ (open 
symbols), is more complicated and will be discussed in Section El 

A baryons. Now let us consider the A baryons. Figure 121 shows the (J, M 2 ) 
plane for all A baryons (filled symbols). The filled circles show the P = +1 baryons 
with smallest masses at a given J (their PDG names are indicated). It should be 
noted that the PDG name for the A(2000) baryon with J = 5/2 is likely to be 
inappropriate since most recent experimental studies have indicated that its mass 
is either 1724 ± 61 MeV [12] or 1752 ± 32 MeV [13]. For the linear fit, we use the 
former mass, which is also quoted in [11]. 

The linear least-squares regression for the lowest-mass A baryons with P = 
+1 is shown in Fig. 12 by the solid line. The ^/ndf = 6.1/4 supports the linear 
fit, despite very small experimental uncertainties on the measured masses and the 
significant range in J. The quality of the linear regression is impressive: The dashed 
lines indicate a 95% confidence-level interval for the linear regression, which is only 
170 MeV wide even for J = 21/2. Both multiple r-squared and adjusted r-squared 
are 0.9999 and the p- value is 6.2- 10~ 9 . Furthermore, a possible non-linear trend was 
checked by using a second-order polynomial function. Such a fit has /ndf = 5.4/3, 
while the parameter for the quadratic term was consistent with zero (0.014 ±0.017). 
Thus, one concludes again that there is no evidence for a non-linear behaviour. 

The linear regression has only one significant residual at J = 5/2 where the 
measured mass squared is 100 MeV 2 above the upper 95% confidence-level line. 
However, as was noted above, the existence of the A(2000) baryon is not well settled 
(two PDG stars), and its mass needs to be determined more accurately. It should also 
be noted that only three baryons shown in Fig. El A(1232), A(1950) and A(2420), 
are fairly well studied. 

Using the principles outline above, it is impossible to find a sufficient number 
of A baryons with P = — 1 for the linear regression fit. This will be discussed in 
Sectional 
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A baryons. Figure El shows the (J, M 2 ) plane for all PDG A baryons (filled 
symbols). The filled circles show the baryons with the smallest masses and with P = 
— 1. Again, such baryons are well described by the linear regression fit (x 2 /ndf = 
0.2/1). 

Figure H3 also shows the A baryons with smallest masses (their names are not indi- 
cated) andP = +1. The latter baryons (A(1115)(P01), A(1890)(P03), A(1820)(P05) 
and A(2350)(PT09)) cannot be used for the fit since: 1) A(1115)(P01) decays weakly; 
2) the position with J = 3/2 is already filled with the A(1520)(P03) baryon which 
has the opposite parity. We examine this further in Sect. El 

Other baryons. For Q and H baryons, the existing experimental data are 
insufficient for the construction of the trajectories with more than two data points. 
The £ baryons will be discussed below. 

3 Discussion 

It is rather convincing now that the method proposed at the beginning of Sect. [21 
indeed works rather well: it allows the construction of three perfectly linear baryonic 
trajectories with more than two data points. The fit parameters are summarised in 
Table □ 

Now, let us discuss what would happen if: 1) weakly decaying baryons will be 
included in the fit as well; 2) one considers baryonic trajectories for a certain parity 
even when there are baryons with the opposite parity with lower masses. In this 
context, we will also discuss hadrons with different space parity which have been 
omitted from the consideration in Section |21 We will not analyse the minimum- mass 
requirement itself since, without it, the vast majority of baryonic trajectories will 
be non-linear and identical to those studied elsewhere [5,7-9]. 

First of all, let us include the lowest-mass weakly decaying A state shown by the 
open symbol at J = 1/2 in Figure El We will exclude the A(1820)(P05) state at 
J = 3/2, which is also indicated with the open symbol, since this hadron does not 
satisfy to the minimum-mass criteria. The linear fit of the remaining three baryons, 
A(1115)(P01), A(1890)(P03) and A(2350)(PD9), cannot be considered as a perfect 
one since x 2 / n df = 3.2/1. Thus, the weakly decaying A(1115)(P01) state violates 
the linearity. 

There is a similar situation with the S baryons. Figure 0] shows the £ baryons 
with P = +1. There are only two points which can be used to construct the 
trajectory. If the lowest-mass weakly decaying E(P11) at J = 1/2 is included, the 
data points cannot be described by the linear fit. 

The open circles in Fig. 0] show two £ baryons whose angular momenta are 
unknown, but their masses are reasonably well determined. The values of the angular 
momenta were hypothesized assuming that such baryons should by close to the line 
determined by £(1385) (P13) and £(2030) (P17) baryons. At the same time, these 
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new baryons should be above the A-trajectory which bounds the entire baryonic 
spectra at low masses. 

Now let us assume that the minimum-mass requirement for a given parity is 
not essential in cases when there exist baryons with the opposite parity but with 
smaller masses for a fixed J. First, let us consider the N baryons with P = — 1 
shown in Fig. QJopen symbols). For J = 3/2 and 11/2, one can use iV(1520)(-D13) 
and iV(2600)(/lll) baryons without any ambiguity. For J = 5/2 and J = 9/2, 
iV(1675)(Z}15) and iV(2250)(G19) baryons with P = — 1 overlap in masses with 
the N baryons of the opposite parity, i.e. iV(1680)(F15) and N(2220)(H 19). If 
one ignores the minimum- mass requirement, then all baryons indicated by the open 
symbols should be considered for the regression. Obviously, the linear fit will fail in 
this case. However, a linear trajectory for the P = — 1 baryons may still exist if the 
iV(1675)(-D15) and A^(2250)(G19) states are removed from the consideration on the 
basis of the minimum-mass requirement. 

Now let us consider strongly decaying A with P = +1 shown in Figure El It is 
evident that the data points shown with the open symbols cannot be described by 
a linear function. Thus, it is essential to exclude A(1890)(P03) at J = 3/2 from 
the consideration. This can indeed be done taking into account the minimum mass 
requirement and noting the presence of the low-mass A(1520)(-D03) state. 

There is another example: A baryons with P = — 1 (see Fig. |3J). Such baryons 
cannot lie on the same line since the linear fit is characterised by x 2 /ndf = 11/3 
and a wide 95% confidence interval. This can be explained as before: positions with 
J = 3/2, 5/2 and 9/2 have already been filled with the P = +1 baryons with lower 
masses (shown in Fig. El with small filled symbols). Therefore, there are only two 
baryons left with J = 1/2 and 11/2, which is insufficient for the linear regression 
analysis. 

4 Towards extraction of the diquark mass 

If our hypothesis is correct, then the present experimental data can be used for 
the construction only three linear trajectories with more than two data points. All 
other baryons lying above such trajectories on the (J, M 2 ) plane have additional 
contributions to their masses from non-pure orbital rotations. Perhaps one can use 
the term "intrinsic noise" [4] for such states: all such hadrons may have a non-linear 
relationship between M 2 and J [5,7-9]. 

The linear trajectories proposed above can be used for the validation of a rela- 
tivistic model in which a hadron can be treated as a rotating flux tube (or a string) 
with a quark and a diquark at the ends. Such a string can be characterised by a 
constant tension T = a /2tt. For small (di) quark masses, an approximate solution 
of this model is given by the Selem-Wilczek (SW) expression [4]: 




3/2 



= m 



3/2 



1 
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where m x and m 2 are the masses of diquark and quark connected by a relativistic 
string. The fitting function directly follows from Eq. (J3J): 

M 2 = a J + \ ■ (crj) 1/4 + \^L. (4) 

3 v 9 v^J 

This equation resembles the Chew-Frautschi relationship between M 2 and J in the 
limit of small masses. 

In order to test the model above, we will use the N baryons shown in Figure ^ 
The fit result using the SW function Eq. (0J) is shown in Fig. |U] (solid line), while 
the two dashed lines illustrate the fit sensitivity to fi (in this case we keep the 
slope parameter a to be the same as for the nominal fit using Eq. (0J)). The dotted 
line shows the linear fit as in Figure The fit with the function given in Eq. (@J) 
is strikingly good. The fit parameters are fully constrained by the proton mass 
at J = 1/2, since this is exactly the region with a highest sensitivity to fi. The 
parameter a is 0.908 GeV 2 , which is close to the slope value obtained from the linear 
fit shown in Figure ^ The extracted value of the mass parameter is fi = 0.323 GeV. 

Similar fits have been performed using the A and A baryons shown in Figures El 
and H3 For such baryons, the fit with Eq. (J3J) did not converge: the parameter /i 
had a fitting error larger than its value and the fit had a significant sensitivity to 
initial fit values. The reason for this is rather simple: As seen from Fig. |UJ the last 
non-linear term in Eq. (0J) can only be constrained by the region J < 1. However, 
A and A trajectories do not have a hadron at such a small J. Even although such 
baryons cannot constrain the mass parameter /i, their linear trajectories are still 
consistent with a fairly linear behaviour of Eq. (@J) at large J. 

It should be noted that the solution given in Eq. (jSJ) was obtained in the limit 
of negligible quark masses (/i — > 0) and for large J. However, the obtained values 
of n and the range of J used to fit the data may lead to a worry that the above 
assumptions are not appropriate and the terms beyond 0(fx 3 ^ 2 ) are necessary to 
consider in Eq. (jSJ). Therefore, we have made an attempt to solve the equations 
of the diquark string model analytically by taking into account high-order terms 
neglected in the solution [4]. We have obtained a relationship between the mass and 
J using a full set of 0(/i 3//2 ) terms and, in addition, some terms of order 0((i 5 / 2 ) 
(see Appendix B). It should be stressed that a complete 0(/i 5,/2 ) calculation requires 
a solution of the appropriate quintic equation. 

Figure [7| shows the fit using the 0(/i 5 ' 2 ) solution given in Eq. (jB-10|) of Ap- 
pendix B. The fit, shown with the dashed line, was performed using three free 
parameters, mi, m 2 and a. It was assumed that the mass of the diquark (mi) is 
larger than that of the quark (m 2 ) during the fit. The quality of the fit is fair. It 
should be noted that if only two parameters are used for the fit, i.e. [i and a as 
in the SW case, x 2 /ndf is smaller (= 2.2/2). In is interesting to observe that M 2 
increases when J decreases. In fact, this is an artifact of truncation of the series in 
Eq. (jB-ll|) : high-order terms proportional to 0(fx 7 ^ 2 ) and 0(fi 9 ^ 2 ) have negatived 
contributions and thus turn to reduce the increase of M 2 at small J. 
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It is rather clear that in order to obtain a reliable model prediction when the 
masses are not too small, it is essential to find a numerical solution of Eq. (jB-5|) . 
The result of our numerical calculation 2 used in a ^-minimisation procedure for the 
iV-baryonic trajectory is shown in Figure [7| (the solid line). The fit was found to be 
excellent. 

Table [21 summarizes the fit parameters obtained using different approximations. 
The numerical solution leads to very similar masses as for the analytical 0(/i 5//2 ) 
calculation. Moreover, the parameter \i ~ 0.314 GeV calculated from the extracted 
mi and m 2 is very similar to that from the SW fit shown in Figure |H1 As before, 
the fit can only be constrained by the J < 1 region, thus the iV-baryon trajectory 
is the most useful for the extraction of the mass parameters and for the validation 
of this model. 

A few words about the precision on the extracted masses are necessary. Table El 
shows the fit values with the necessary numerical precision to reproduce the x 2 /ndf 
and thus the proton mass, which has a very small experimental uncertainty. The 
uncertainties on the extracted parameters from the fits are negligible, therefore, they 
are not quoted (see Appendix A for details). We did not estimate the exact range 
of the fit parameters which give an acceptable fit (i.e. with x 2 jndj < 1) since this 
will require a significant computational time. However, from several tests we did, 
we have concluded that the parameters are fairly close to ni\ ~ 0.27 ± 0.01 GeV, 
m 2 ~ 0.11 ± 0.01 GeV and a ~ 0.92 ± 0.1 GeV 2 (the so-called Fit I), with a very 
strong and complicated correlation between the values. 

In addition to the solution given above, we have found another solution with 
mi = 0.228 GeV, m 2 = 0.179 GeV and a = 0.920 GeV 2 (Fit II, see Table C-l of 
Appendix C), which also gives an acceptable fit, x 1 /ndf = 0.9/1. In this case, both 
masses are rather similar, almost in the spirit of the expectations discussed in [4] . 
This solution corresponds to a different radius of the quark-diquark system (see a 
discussion in Appendix C). We would like to note that the obtained parameters 
may not be the only solutions which lead an acceptable fit with x 2 /ndf < 1, thus 
the model uncertainties are unknown yet. To find possible alternative solutions will 
require a significant computational time and, therefore, this is outside of the scope 
of this paper. The model uncertainties for the iV-baryon trajectory are further 
discussed in Appendix C. 

It was also verified that if one sets mi = m 2 for the fit using the numerical 
solution, then the fit cannot be considered as a good one, x 2 / n df — 5/2. Obviously, 
the assumption mi = m 2 is more appropriate for mesons, and it is remarkable that 
the iV-baryon trajectory does not support it. 

2 We used the function DZERO from the CERNLIB FORTRAN library [14] to find a zero of a 
real- valued function for solving Eq. I|B-5(I . Then such a numerical solution was used in the fitting 
function using the MINUIT program, see Appendix A. 
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5 Mesons 



It would be interesting to see the applicability of Eq. (J2J) and the assumptions 
discussed in Sect. El to the mesonic sector. We will mostly be interested in the 
linearity of the trajectory which bounds all known mesons on the (J, M 2 ) plane at 
low masses. This is very similar to the A baryons which bound the entire baryonic 
sector at low masses by the straight trajectory shown in Fig. |21 

Figure |H1 shows the (J, M 2 ) plane for all PDG mesons (filled symbols). As before, 
the filled circles show all mesons with lowest masses (for J > 1) which were used 
in the linear fit. It can be seen that such mesons approximately lie on a straight 
line, and this can scarcely be a coincidence. The fit has a significant x 2 / n df due to 
very small experimental uncertainties on the measured meson masses, thus one can 
judge about the linearity only with a certain caution. There is only one significant 
residual at J = 5; however, it is not improbable that the A(2210) state may have 
J = 5 and thus it falls exactly on the linear fit. All other states indicated with the 
open circles have unknown J and their locations were hypothesized. 

While the slope for the meson sector is somewhat lower, it should be noted 
that the upper 95% confidence-level line (po = — 0.503, pi = 1.099) approximately 
coincides with that for the baryonic sector. 

It is interesting to note that there is a clear periodicity in space and charge 
conjugation parity indicated as (P, C) for lowest-mass mesons which bound the 
meson spectrum on the (J, M 2 ) plane (see Fig. |HJ). This was already noted for the 
trajectory of vector mesons in [3]. Based on this observation, we can predict that 
the C-parity of A (2750) is likely to be negative. 

6 Conclusion 

In this paper, using the most recent experimental data, empirical rules have been 
proposed which are sufficient to reveal a strict linear relationship between the total 
angular momentum and squared-mass of baryons. Using a least-squares regression, 
we did not find evidence for non-linearity of such trajectories even when the most 
stringent statistical tests based on the \ 2 an d 95% confidence-level intervals were 
used. In one case, this conclusion was made after the analysis of a number of baryons 
as large as six. This observation may provide the basis for a new systematization of 
hadrons and certainly require thoughtful theoretical investigation. 

The linear A-baryon trajectory can be well described by a semi-classical approach 
based on a spinning relativistic string with diquark and quark at the ends. We have 
shown this baryonic trajectory is \ 2 consistent with the exact numerical solution of 
this model (up to nine digits of the extracted model parameters). In this paper we 
have determined the (di)quark masses from the fit of the A-baryon trajectory using 
our solution. 

In principle, this semi-classical model is qualitatively consistent with other linear 
baryonic trajectories. However, in this paper, we did not attempt to verify this 
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using our numerical calculations. Also, we did not analyse mesons as carefully as 
baryons. However, it is interesting to observe that the meson spectrum on the 
(J, M 2 ) plane is restricted at small masses by a linear trajectory formed by mesons 
with S = C = B = which have periodicity in space and charge conjugation parity. 
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APPENDIX A 



The linear regression was performed using the weights w = 1/<t 2 , where a is the 
experimental uncertainty on the mass squared of a hadron. In case of asymmetric 
experimental errors, we take the average of the upper and the lower experimental 
uncertainty 

The least-square linear regression was carried out using the R program [15]. The 
R-package was also used to estimate confidence-level intervals on the linear regres- 
sion. This package gives larger uncertainties on the linear-regression parameters 
than the MINUIT program from the CERNLIB FORTRAN library [14]. The MI- 
NUIT parameter errors give information on the uncertainty in the best fit values 
and are not meaningful when points do not have (or have very small) experimental 
uncertainties. 

In contrast to the standard linear least-square regression analysis, the MINUIT 
program was used for validation of the relativistic string model discussed in Sec- 
tion This simplifies the fit procedure in case of non-linear functions, especially 
when their analytic form is unknown (as in case of the numerical solution of the 
string model discussed in Section HJ). The extracted fit parameters with the neces- 
sary numerical precision to reproduce the proton mass are shown in Table El The 
MINUIT parameter errors are not shown, since they are smaller than the last digit. 



APPENDIX B 

Here we will derive the functions used to fit the iV-baryon trajectory shown in Fig. [71 
First, recall the main set of equations (8)-(9) of the string model [4]: 

T 



E = ^ ( m «7n + — arcsin(/5 n )j, (B-l) 

n=l 



n=l 



T = cjm„/5 n 7 n , n = l,2. (B-3) 

Here, E is the energy of the quark-diquark system rotating with an angular momen- 
tum J, m n denotes the mass of a diquark (n — 1) and quark (n = 2), u is angular 
velocity, T is a string tension and f3 n is a linear velocity. The factor 7 n of a given 
(di) quark is defined as 

1 

where r n is a distance from the center of rotation. 

Our goal is to find a relationship between the energy E and J. Such a relationship 
has to be expressed in terms of a = 2ixT and m n . Excluding terms with arcsin(/3 n ) 
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from Eq. (|B-1|) and using the representation 



• (a \ n 1 • fo@ n 
arcsm (p n ) = arcsm 2 — 



one gets from Eq. (jB-l|) - (jB-3 

E : 



mi mo 
2uJ + — + — . 



7i 



72 



2T\uj) 
1 



arcsm 



n=l 



2Pn 
In 



J- /. 



n = 1,2. 



(B-4) 
(B-5) 
(B-6) 



We use Eqs. flB-4|) - ()B-6|) to fit the A^-baryon trajectory and to extract the parameters 
m n and o from the data. At fixed J, m n and a, the equation (jB-5|) can be solved 
relative to u either numerically, or to some approximation, using a series expansion. 
Here, x n = m n uj/T can be used as the expansion variable in the limit of small 
(di)quark masses. 
In terms of x„, 



1 

In 

Pn 
In 



*2x n 



x n + a/4 + x\ 



1/2 - T l/2 _ i_ T 3/2 , * 5/2 , a(JT/2\ 
■^n n ? n ~ nfi n ~ ^V^n J 



2 5 * 



3/2 



Using the expansion 



T V2 _ Ji T 3/2 

n 2 2 



9 
25 



x^ /2 + £>(x 



7/2 \ 



(B-7) 
(B-8) 



arcsin ( 2 



7n -3 V 7n 



5 V7n 



Eq. (|B-5|) becomes 



2TL 



7T 



2/#l\3 2//3 2 \3 6 /A 



3 V71 



3 V72 



5 V71 



6 ffc 



5 \7 2 



+ 0( 



PnV 
In 



(B-9) 



Thus, substituting Eq. (jFTTj) and Eq. (jFT8jl in Eq. (jB^4jl and Eq. (jFT9jl and keeping 
only terms up to order m^ 2 , one obtains: 



E 
J 



7T 



T 



CO 



T) 20 VT/ 



+ 0(/i 7/2 ), 



(B-10) 



JL (1\ 2 _ $1 C£\ ~ 1/2 + jLm1(^\ 1/2 



2T\u 



3T VT 



20 T VT 



C(^ 7/2 ) (B-ll) 



.,, 3/2 3/2 . 3/2 , 5/2 5/2 . 5/2 
with /i 1 = m{ + m 2 and /j 2 = + m 2 . 
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Let us denote V = \/ujJT and rewrite Eq. (|B-11|) as a quartic equation 

V A + aV 3 - b = 0, (B-12) 

with 



3/2 



3aJ 

5/2 



a J 



where the term proportional to // 2 was neglected. To solve the quartic equation 
(IB- 12)1 we use the Ferrari method [16]. The resolvent cubic of Eq. ()B-12|) 



fT + 46(7 + ba 2 = 0, 
transforms Eq. ()B-12)) to a quadratic equation: 

v 2 + 7± y + 5± = 

with 



7± = -[a±Va 2 + 4C/], <J ± = -[U ± y/U 2 + 4b], 
where a real root of Eq. ()B-13)) (the Cardano formula) is 



U =- 
6 



126V810 4 + 7686 - 1086a 2 



1/3 i 

~ 6 



126V810 4 + 7686 + 1086a 2 



1/3 



(B-13) 

(B-14) 
(B-15) 

. (B-16) 



Thus, a real and positive solution of Eq. ()B-14)) and therefore, of Eq. (|B-12)1 . is 

a 2 +iuy' 2 - 



T 



+ 



(a 2 + 4(7) 1 /2 - a 



+ 



{U 2 + 4b) 1 / 2 - U 

(B-17) 

Substituting the above expression in Eq. (|B-10)I . one obtains a desired relation- 
ship between the baryonic mass M = E and variables J, a and m n . 



APPENDIX C 

Using a graphical representation of different parts of Eqs. ()B-4)1 and ()B-5jl . it is 
easy to demonstrate that the string model described by the system of equations 
Eqs. ()B-4)) - (|B-6)) allows several solutions for the (di) quark masses. Such solutions 
have trajectories with very small (but different) non-linear behaviour. 
Let us denote by Z n the following variables 

Z n = ^, n=l,2 (C-l) 

In 

and introduce a function 

2 

F(Z X , Z 2 ) = J2 (Zn - arcsin(Z n )) (C-2) 

n=l 
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which represents the last term in Eq. (jB-5|) . The variables Z n have values in the 
interval [0, 1], and are related to the parameters m n , uj and T through the equations 
Eq. (jB-6|) and Eq. (jB-8|) . The function F(Zi, Z 2 ) takes values in the range [2 — ir, 0]. 
As shown in Fig. IC-lf a). F(Zi,Z2) is symmetric with respect to the median line 
Zi = Z2 and a monotonic function of the both arguments. The main feature of 
Eq. (jB-5|) is that the equation F(Z 1; Z 2 ) = J 7 , where J 7 is a constant, defines an isoline 
on the surface and at fixed J all the points {Zi,Z 2 } on this isoline correspond to 
the same angular velocity uj. In Fig. IC-lf a). the isolines are shown as lines between 
colored bands. 



F(Z, ,Z 2 ) 




Figure C-l: F(Z 1; Z 2 ) as a function of Z x and Z 2 (left) and Z n shown as a function 
of the (di)quark velocity (right). 

Let us consider the mass evolution equation (|B-4|) starting from a very high J. 
As follows from (jB-5|) and the definition of x n , a very high J can be obtained only 
at very small uj or x n . Then Eqs. (B-7) and (B-8) tells that 7 n 3> 1, (3 n « 1, Z n « 
and F(Zi, Z 2 ) « 0. Thus, at such J, /3„ and uj, the quark and di- quark are at large 
distances from the center of rotation. 

Now, if J decreases and uj, consequently, increases, the phase point (Zi,Z 2 ) on 
the surface F(Z 1; Z 2 ) is moving along some path in the direction of (Zi = 1, Z 2 = 1). 
A fit of the A^-baryon trajectory fixes a particular path passing at J = 1/2 through 
a point (Zi^, Z 2 ^) on the surface given by F(Z 1; Z 2 ). If Z n function is plotted as 
a function of the (di)quark velocity /3 n , as shown in Fig. IC-l( b). one finds that the 
point is moving from (3 « 1 to a smaller (5. The maximum Z n = 1 is reached at 
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Numerical solution 


Numerical solution 


limit L\ — L'i — 1, 


Parameters 


TP if T 

r it 1 


TP if TT 

r it 11 




-1,-2 / /y. r$ -p 

X l na J 


u.y/ 1 


n n/1 
u.y/ 1 




777 -i (tP V 


D 971 71 3 


228808 


18393 


777 o 


0.114134 


179558 


18393 


a, GeV 2 


0.923374 


0.920000 


1.12428 


J 


1/2 

/ 


1/2 


1/2 

/ 


u, GeV 


0.643476 


0.628674 


0.679332 




0.935557 


0.974606 


1 




0.975607 


0.998309 


1 


F^Z^) 


-0.648139 


-0.884682 


2-7T 




1.189721 


0.982401 


1/V2 


^2 


0.499746 


0.770943 


1/V2 




0.822549 


0.782280 


1/V2 




0.780872 


0.727369 


1/V2 


rx, GeV" 1 


1.278291 


1.244335 


1.045218 


r 2 , GeV" 1 


1.213521 


1.156988 


1.045218 



Table C-l: The model parameters extracted from Eqs. (jB-6l) - (IB-8 



A) = VV2. 

Table C-l gives the model parameters 3 extracted from Eqs. ()B-6j) - (jB-8)) using 
m n and cr found from Fit I and Fit II. First of all, let us note that all (3 n < 1/ y/2 
(or higher values of uS) correspond to a non-physical J < 1/2 and (5 n > l/v2 holds 
for the both fits. 

The baryon masses predicted by the string model are degenerated with respect 
to uj and m n . Indeed, F(Z 1; Z 2 ) is negative and, therefore, there is a range of ui for 
which one obtains J\ = J 2 at different u)\ and uj 2 . If (B-5) is equated at different u, 
then 



-i[2 7 r + ^ 1 ] = ^[27r + ^ 2 ] 



or 



-F 2 = (^) 2 (2tt + JF x ) _2tt. 



(C-3) 



(C-4) 



Here T\ and Ti denotes F(Zi, Z 2 ) at two phase-space points (Zi, Z 2 )i )2 . The condi- 
tion 

2 - 7t < JS < (C-5) 
puts the following restriction on aj 2 (at a fixed lo\): 



2 + 7T 
2^ + ^! 



< 



Co> 2 



< 



2tt 



2n + F 1 



(C-6) 



3 We have reduced the numerical precision of the quoted fit values for representative purposes. 
The exact fit values can be found from Table 2. 
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The condition (|C-6|) could mean that if one finds a set of parameters m 2 , a} 1 
with the evolution along a path {u, /3 1; /^li, there should exist another solution 
{mi, m 2 , cr} 2 with a different path {u, Pi, /3 2 } 2 . The second found solution for the 
iV-baryon trajectory may confirm such interpretation. Thus, we conclude that the 
model should have an intrinsic uncertainty for the extracted masses depending on 
the Zi and Z 2 values. 
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Particle family 


po, GeV 2 


pi, GeV 2 


X 2 /ndf 


iV baryons, P = +1 
A baryons, P = +1 
A baryons, P = — 1 
E baryons, P = +1 
mesons (lowest mass) 


0.393 ±0.003 
-0.156 ±0.013 
0.746 ± 0.007 
0.277 ±0.024 
-0.454 ±0.026 


0.974 ± 0.006 
1.116 ±0.006 
1.042 ± 0.004 
1.098 ±0.016 
1.052 ±0.019 


0.5/2 
6.1/4 
0.2/1 

220/5 



Table 1: Fit parameters obtained using the linear parameterisation M 2 = po + p\J 
for several hadronic families with the lowest masses on the (J,M 2 ) plane. The 
X 2 ' /ndf for the E states is not given since there are only two points for the linear fit 
(see the text). 



Parameters 


SW fit 


0(^ 2 ) 
solution 


G(/i 5 / 2 ) 
solution 


numerical 
solution 


a, GeV 2 
/i, GeV 
mi, GeV 
m 2 , GeV 
X 2 /ndf 


0.908091023 
0.323158469 

1.5/2 


0.889953024 

0.302678728 
0.109731744 
2.2/1 


0.902776942 

0.280596424 
0.109180448 
1.8/1 


0.923374256 

0.271713077 
0.114133666 
0.9/1 


M p , GeV 
from Eq. ()B-4J) 


0.938272037 


0.938272031 


0.938272034 


0.938272048 


M p , GeV [10] 


0.938272029 ± 0.000000080 



Table 2: Fit parameters obtained using different solutions of the relativistic diquark 
model (see the text). The A^-baryon trajectory was used for the fit. The parameter 
uncertainties are not shown (see Appendix A). The large number of digits, which is 
necessary to reproduce the quoted x 2 / n df, is due to a large experimental precision 
on the A^-baryon masses. The proton mass was calculated from Eq. (jB-4|) using 
different approximations and compared with the world average value. 
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N baryons 



M (GeV) 




J 

Figure 1: The (J, M 2 ) plane for iV-baryons. The linear fit was performed using 
the P = +1 baryons with smallest masses for a given J. The solid line shows the 
weighted-least squares regression using a linear model and the dashed lines show 
the corresponding 95% confidence- level interval. The open symbols show lowest- 
mass baryons with P = -1: JV(1520)(D13) (with J = 3/2), iV(1675)(D15) (with 
J = 5/2), JV(2190)(G17) (with J = 7/2), iV(2250)(G19) (with J = 9/2) and 
JV(2600)(nil) (with J = 13/2). The JV(1675)(D15) state overlaps in mass with 
the iV(1680)(-F15) baryon with positive parity. 
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A baryons 



M (GeV) 




J 



Figure 2: The (J, M 2 ) plane for A baryons with P = +1. The filled circles show 
lowest-mass A baryons used for the linear weighted-least squares regression. The 
solid line shows the fit, while the dashed line indicates a 95% confidence-level inter- 
val. Small filled squares show all other A baryons. 
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8 
7 
6 
5 
4 
3 
2 
1 




M (GeV) 



A baryons 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



2.742 



2.545 



• Used for the fit (P=-1) 
O Lowest mass A with P=+1 

Linear fit 

95% confidence level 




M =pO+p1*J 
X 2 /ndf = 0.2/1 
p0 = 0.746 ±0.007 
p1 = 1.042 ±0.004 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

J 



Figure 3: The (J, M 2 ) plane for A baryons with P = — 1 (filled symbols) and P = +1 
(open symbols). The latter hadrons correspond to A(1115)(P01), A(1820)(F05), 
A(1890)(P03) and A(2350)(#09). The A(1820)(F05) baryon with P = +1 is on top 
of A(1830)(-D05) with P = —1. Small filled squares show all other PDG A baryons. 



20 



£ baryons 




i i i i i i i i ■ i i i i i ■ i i i i i ■ i i i i i ■ i i i i i ■ i i i i i ■ i i i i i ■ i i i i i ■ i i i i i ■ i i i i i i i 
123456789 10 11 



J 

Figure 4: The (J, M 2 ) plane for E baryons. The filled circles show the S baryons 
with P = +1. The filled triangle shows the lowest-mass weakly-decaying £ with 
P = +1, while the small filled squares show all other PDG E baryons. The open 
circles indicate possible locations of several high-mass PDG baryons with unknown 
J. 
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10c 



M (GeV) 



A baryons 



i i i i l i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i, 'i i i i i i i i i i i i i 



• Used for the fit (P=-1) 
Linear fit 

95% confidence level 




Figure 5: As Fig. |2l but the fit was performed for the lowest-mass A baryons with 
P = -l. 
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Figure 6: The fit results for leading and daughter nucleons with positive parity. 
Only baryons with the minimum mass requirement for a given J were used. The 
solid line indicated the fit by using Eq. (J1J (the so-called SW fit) with two free 
parameters, \x and a. The parameter uncertainties are negligible (not shown). In 
order to illustrate the fit sensitivity to fi, the dashed lines show the same fit function 
but with different \i. The dotted line shows the linear fit shown in Fig. ^ 
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t 

\ 
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fit using 0(n 5/2 ) solution 








t 
% 
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% 


numerical solution 








1 

I 
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« 

t 






f\ ^(1680)(F15) 




t 

\ 

— * 

X 




P(P11) 








v^-— — T» 


0(n 5/2 ) solution: 




numerical solution: 








X 2 /ndf = 1.8/1 




X 2 /ndf = 0.9/1 








= 0.28060 




m 1 = 0.27171 








m 2 = 0.10918 




m 2 = 0.11413 








a =0.90278 




a =0.92337 




1 


I i i 


i I 


I 





l"1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

10 2 10" 1 1 10 



log J 

Figure 7: The dotted line shows the SW fit as in Figure |Hl The dashed line shows 
the fit using the solution of the diquark string model calculated up to 0(/i 5 / 2 ) terms, 
see Eq. (IB-10|) of Appendix B. The fit was performed using two mass parameters, 
rrii and m 2 (wii > m 2 ) and a string tension a. The solid line shows the fit using the 
exact numerical solution (see the text). The parameter uncertainties are negligible. 
Table El gives more exact numbers from the fits. 
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Mesons 

M(GeV) | 




123456789 10 11 

J 



Figure 8: The (J, M 2 ) plane for mesons. The filled circles show the mesons with 
smallest masses and J > 1 used for a weighted-least squares linear regression (the 
solid line). The dashed lines indicate a 95% confidence-level interval for the linear 
regression. For each meson used in the fit, space parity (P) and charge conjugation 
parity (C) are indicated in parentheses. Small filled squares show all other PDG 
mesons and the open circles indicate guessed location of several high-mass mesons 
with unknown J. 
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